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Abstract 

The melting behavior of an Ag microwire mesh induced by Joule heating was numerically investigated and 
compared with that of the corresponding Ag nanowire mesh with the same structure but different geometrical 
and physical properties of the wire itself. According to the relationship of melting current and melting voltage 
during the melting process, a similar repetitive zigzag pattern in melting behavior was discovered in both meshes. 
On this basis, a dimensionless parameter defined as figure of merit was proposed to characterize the current-carrying 
ability of the mesh. The consistent feature of figure of merit in both meshes indicates that the melting behavior of the 
Ag nanowire mesh can be predicted from the present results of the corresponding Ag microwire mesh with the same 
structure but made from a different wire (e.g., different size, different material) through simple conversion. The present 
findings can provide fundamental insight into the reliability analysis on the metallic nanowire mesh-based transparent 
conductive electrode. 

Keywords: Consistent feature; Figure of merit; Joule heating; Melting behavior; Ag microwire/nanowire mesh; 
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Background 

To meet the requirement of next-generation flexible 
optoelectronics for both information (e.g., display, elec- 
tronic reader, touch screen) and energy (e.g., solar cell, 
window glass), there is growing interest to develop 
transparent conductive electrodes (TCEs) possessing 
high optical transmission, good electrical conductivity, 
and excellent flexibility [1,2]. However, the present com- 
mon commercial TCE material, i.e., indium tin oxide 
(ITO), suffers from several major limitations [3-5], such 
as high cost due to the shortage of indium and poor 
mechanical stability due to the brittleness. Therefore, 
it is highly desirable to find a promising alternative 
which can be used in the forthcoming TCEs [6]. Re- 
cently, the network of various nanostructured materials 
(e.g., carbon nanotube [7,8], graphene [9-11], metallic 
nanowire [12-20] /nanotrough [21] /honeycomb [22], 
and the combinations of the above [3,23-25]) has shown 
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great potential for the application in optoelectronic devices 
such as solar cells [9,16-18] and touch screens [14,20]. 

Here, our focus is on metallic nanowire mesh (i.e., 
regular nanowire network) because of its ideal character- 
istics of low sheet resistance, high optical transparency, 
and flexible controllability. For example, Kang et al. [16] 
have fabricated a Cu nanowire mesh electrode on a poly- 
ethylene terephthalate (PET) substrate, which shows 
compatible optical transmittance in the visible wave- 
length range with commercial ITO-coated PET and offers 
lower sheet resistance than ITO. Moreover, the short- 
circuit current and power conversion efficiency of a solar 
cell with this type of Cu nanowire mesh electrode are 
comparable to those of the same device using an ITO 
electrode. 

However, realizing the potential benefits of such me- 
tallic nanowire mesh in practical optoelectronic devices 
remains a great challenge because of the lack of reliabil- 
ity analysis. It is known that the pathway of current in a 
metallic nanowire mesh remains in the nanowire itself, 
instead of uniform distribution throughout the whole 
ITO film. Great reduction in current flow area will cause 
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enormous increase in current density and significant rise 
in temperature due to Joule heating. Therefore, it is be- 
lieved that the melting induced by Joule heating is a po- 
tential threat to the degradation of the metallic nanowire 
mesh-based TCE, which may cause deterioration of the 
corresponding optoelectronic devices. In a pioneering 
experimental report, Khaligh and Goldthorpe [26] have 
indicated that at a constant current density, a random 
Ag nanowire network fails after a certain period. More- 
over, the network with higher sheet resistance carrying 
greater current density will fail more easily because of 
Joule heating. Hereafter, a numerical method has been 
developed [27] by the present authors to clarify the melt- 
ing behavior of metallic nanowire mesh due to Joule 
heating. Using this technique, a repetitive zigzag pattern 
in the relationship of melting current and melting volt- 
age triggering the melting of the mesh was discovered. 
It indicates that in real working conditions, a metallic 
nanowire mesh supplied with current source may ex- 
perience repetitive unstable (where several wires are 
melted simultaneously at a constant current/voltage) and 
stable (where an increment of current/voltage is neces- 
sary for melting progression) melting behavior until the 
mesh is open. However, some of these predicted intrinsic 
features in the melting of the metallic nanowire mesh 
would not be detectable because of the difficulty in sam- 
ple preparation and experimental measurement. 

To overcome the above weakness, the relatively easy- 
to-prepare microwire mesh comes into the sight. One 
might expect the melting behavior of microwire and 
nanowire meshes to be similar by assuming that the cur- 
rents would just scale up. However, metallic nanowire in 
general displays different properties from microwire be- 
cause of significant size effect. For example, with decreas- 
ing dimension, melting point and thermal conductivity 
decrease while electrical resistivity increases. Such differ- 
ences make it difficult to insist on the similarity of the 
melting behavior for microwire and nanowire meshes, even 
if both of which have the same structure under the same 
working conditions. 

Herein, to find the intrinsic relationship of the melting 
behavior between metallic microwire and nanowire 
meshes, the melting behavior of an Ag microwire mesh 
was numerically investigated and compared to that of 
the corresponding Ag nanowire mesh, which has the 
same mesh structure but different geometrical and phys- 
ical properties of the wire itself. A similar zigzag pat- 
tern was observed in the relationship between melting 
current and melting voltage of both meshes. Therefore, a 
dimensionless parameter defined as figure of merit was 
proposed to indicate the current- carrying ability of the 
mesh. The consistent figure of merit during the whole 
melting process of both meshes implies that the melting 
behavior of the nanowire mesh is predictable from that 



of the microwire mesh by simple conversion. The present 
findings provide fundamental insight into the reliability 
analysis on the metallic nanowire mesh hindered by diffi- 
cult sample preparation and experimental measurement, 
which will be helpful to develop ideal metallic nanowire 
mesh-based TCE with considerable reliability. 

Methods 

A previous numerical method [27] was employed to in- 
vestigate the melting behavior of an Ag microwire mesh 
and compared with that of the corresponding nanowire 
mesh which has the same mesh structure (e.g., pitch size, 
segment number, and boundary conditions) but different 
geometrical and physical properties of the wire itself (e.g., 
cross-sectional area, thermal conductivity, electrical resist- 
ivity, and melting point). 

The mesh structure is illustrated in Figure 1. It is a 
regular network with 10 columns and 10 rows, which 
indicates that the mesh size M@N is 10@10. The 
pitch size / is 200 (im, making the mesh area S of 
3.24 x 10 6 |im 2 . A mesh node (z, ;) denoted by integral 
coordinates (0 < i < M - 1, 0 <j < N - 1) is the intersec- 
tion of the (z + l)th column and the (j + l)th row in the 
mesh. A mesh segment is the wire between two adja- 
cent mesh nodes. For simplicity, the segments on the 
left, right, downside, and upside of the mesh node (z, ;) are 
denoted by S^*, S?.*, S?^, and 5^, respectively. Obvi- 
ously, there areMxA[= 100 mesh nodes and M(N- 1)+N 
(M - 1) = 180 mesh segments. 

The electrical boundary conditions are also shown in 
Figure 1. The load current / is input from node (0, 0) 
and is output from node (9, 0) with zero electrical po- 
tential at node (9, 9). Moreover, there is no external in- 
put/output current for all the other nodes. For the 
thermal boundary conditions, the temperature of the 
peripheral nodes (i.e., (z, 0), (0, j), (z, 9), (9, /)) is set at 
room temperature (RT, T 0 = 300 K), while there is no ex- 
ternal input/output heat energy for all the other nodes. 

The geometrical and physical properties of the wires 
are listed in Table 1. Here, A is the cross-sectional area 




Figure 1 Structure of a wire mesh with size of 10@10 and its 
electrical boundary conditions. 
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Table 1 Geometrical and physical properties of the wires 




Ag microwire 


Ag nanowire 


Al nanowire 


Side length, w (urn) 


1.000 


0.1000 


0.1000 


Cross-sectional area, A (x10~ 2 urn 2 ) 


100.0 


1.000 


1.000 


Melting point, 7~ m (x10 3 K) 


1.234 


0.873 [30] (exp.) 


0.736 [31] (num.) 


Thermal conductivity at RT, X (x10" 4 W/um-K) 


4.200 


3.346 [28] (num.) 


1.150 [32] (num.) 


Electrical resistivity at RT, p 0 (x10~ 2 O-um) 


1.590 


1 1 .90 [29] (exp.) 


6.20 [32] (exp.) 


Electrical resistivity at 7~ m , p m (x10~ 2 Q-um) 


7.200 


37.80 


17.72 



calculated from the side length w of the wire with the 
square cross section, T m is the melting point, X is the 
thermal conductivity, and p is the electrical resistivity 
with the subscripts '0' and 'm' representing the value at 
T 0 and r m . Note that p m [=po{l + a(T m - T 0 )}] is calcu- 
lated by using the temperature coefficient of resistivity a. 
Note that the bulk values of Ag were employed for the 
microwire, while size effect was taken into account for 
the nanowire. Because of electron scattering from sur- 
face and grain boundaries [28,29], the thinner the wire 
is, the lower T m , A, and electrical conductivity {-Hp) 
are. 

To clarify the melting behavior of the mesh, the funda- 
mental theoretical analyses [27] on the corresponding 
electrothermal problem is summarized in the following. 
First, as shown in Figure 2a, a horizontal mesh segment 
(i.e., a wire) S^* between node (i - 1, ;) and (z, ;) with an 
electrically and thermally insulated surface was consid- 
ered, where the current flows from node (i - 1, /) to (r, /). 
Based on Ohm's law, the current density / in the mesh 
segment can be calculated as 



Idcp 
p dx 



(i) 



Here, 0 is the electrical potential, and x is the axial 
coordinate in the mesh segment with the direction right- 
ward for horizontal segment and upward for vertical 
segment. Using Fourier's law, the heat flux q in can 
be calculated as 



-X 



dT 

dx ' 



(2) 



where T is temperature. 

By ignoring heat transfer of the mesh to the underlying 
substrate for simplicity, the heat conduction equation 
can be given as 



d 2 T 



0. 



(3) 



Assuming that the temperatures of nodes (i - 1, /) and (z, 
/) are T^.^ (x = 0) and (x = l), temperature distribu- 
tion in the mesh segment S\ t ^ can be obtained by solv- 
ing Equation 3 as 



9 ,' 2 r 2 

2X ]% 



i 



2X 



i 2 ]x 



(4) 
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Note that in the present simulation, p m was used for p 
to approximate real condition neglecting the effect of 
the temperature dependence of electrical resistivity. 

Second, as shown in Figure 2b,c, the current and heat 
energy passing through a mesh node (z, ;) with four adja- 
cent nodes were considered. In Figure 2b, the current 
is assumed to flow rightward in the horizontal direc- 
tion and upward in the vertical direction. According to 
Kirchhoff s current law, we have 



^external H~ -^internal — 0- 



(5) 



Here, /external is the external input/output current at 
node (/, j), and /internal is the sum of internal currents 
flowing through the node (/, /) from its four adjacent 
nodes. By assuming that the current flowing into the 
node is positive and the current flowing out of the node 
is negative, we can obtain 



* internal 



-J sl 



~Js D ' 



~Js v 



(6) 



where the subscript of / denotes the corresponding mesh 
segment. Taking into account a system of linear equations 
for the node (z, j) composed of Equations 1, 5, and 6, the 
current density in any mesh segment can be obtained. 
Similarly, according to the law of the conservation of heat 
energy, we have 



^external 



internal 



o. 



(7) 



Here, Q ex ternai is the external input/output heat energy 
at node (z, ;), and Qintemai is the sum of the internal heat 
energy transferred to node (z, ;) from its four adjacent 
nodes. In Figure 2c, by assuming that the incoming heat 
energy is positive and the outgoing heat energy is nega- 
tive, we have 



^internal 



{ qch -qrR 



(*»/) 



-4s] 



(8) 



Taking into account a system of linear equations for the 
node (z, ;) composed of Equations 2, 7, and 8, the 
temperature at any mesh node can be obtained. Finally, 
by substituting the above obtained current density in 
any mesh segment and temperature at any mesh node 
into Equation 4, the temperature distribution in any 
mesh segment can be monitored. 

A synopsis of the corresponding computational algo- 
rithm [27] is provided as below. Initially, a small value is 
assigned to the input current /. The corresponding max- 
imum temperature in the mesh r max can be identified, 
which rises with the increasing /. By gradually increasing 
/ with increment AI to make r max reach T m , the first 
mesh segment melts and breaks from an arbitrary small 
force occurring in actual operation (e.g., vibration). At 
that time, the input current and the voltage between 
node (0, 0) and node (9, 0) are recorded as melting 



current I m and melting voltage V m . The corresponding 
resistance R m of the mesh can be calculated by dividing 
Kn by I m . It should be noted that AI must be small 
enough so that melting segment can melt one by one as 
far as possible. Subsequently, an ultra-small value is 
assigned to the cross-sectional area of the first melted 
mesh segment in order to approximate zero. The path- 
way of the current and heat in the mesh is therefore 
renewed. By repeating the aforementioned process, the 
current triggering the melting of mesh segment one by 
one can be obtained until the mesh becomes open. 
Therefore, the relationship between I m and V m as well as 
the variation of R m with the number # b of the broken 
mesh segments can be obtained over the entire melting 
process of the mesh. 

Results and discussion 

Melting behavior of the Ag microwire mesh 

As shown in Figure 3a,b, the obtained relationship of 
melting current I m and melting voltage V m as well as the 
variation of mesh resistance R m with the number n h of 
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Figure 3 Comparison of melting process for both meshes. 

(a) The relationship between / m and V mi and (b) the variation of 
with n b . 
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broken mesh segments during the entire melting process 
of the Ag microwire mesh is compared with those of the 
corresponding Ag nanowire mesh, respectively. 

Obviously, a repetitive zigzag pattern is observed in 
the relationship of I m and V m in the Ag microwire mesh, 
which demonstrates the repetition of three different 
trends: increase of both I m and V m , decrease of both I m 
and Vm, and decrease of I m but increase of V m . Such 
pattern in the melting behavior of Ag microwire mesh is 



similar with that of the corresponding Ag nanowire 
mesh [27]. Note that the microwire mesh has higher I m 
but lower V m than the nanowire mesh, because the 
microwire has larger cross-sectional area and lower re- 
sistivity (see Table 1) and therefore lower electrical re- 
sistance (see Figure 3b) than the nanowire. 

During the melting process, the symmetrical mesh 
structures at three special moments for both meshes are 
compared in Figure 4. The difference in the melting 
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Figure 4 Mesh structures at three special moments in the melting process of both meshes, (a) The starting moment, (b) the moment 
with the maximum current (i.e., sudden fall of current), and (c) the ending moment. 
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pathway of both meshes can be attributed to the differ- 
ent AI for monitoring the melting of mesh segment, 
which are 0.1 mA for the Ag micro wire mesh and 
0.1 uA for the Ag nanowire mesh. Note that such differ- 
ence can be removed by employing much smaller AI for 
the Ag microwire mesh at the expense of increasing 
computational cost. 

Moreover, from the present simulation results, it is be- 
lieved that under constant current density (i.e., current- 
controlled current source), electric breakdown of the 
mesh will never happen as long as the load current / 
does not reach the maximum value of I m (i.e., 7 mC ) even 
if several mesh segments melt. This point is quite differ- 
ent from the reported electrical failure of a random Ag 
nanowire network [26] under constant current density 
after a certain current stressing period. Such difference 
between experiments and present simulations also implies 
that the electrical failure in real Ag nanowire mesh should 
be the synergy of Joule heating and some other possible 
causes, such as corrosion by sulfur, atomic diffusion in the 
nanowire itself, and Rayleigh instability [26] . 

Proposal of figure of merit Z 

To explore the intrinsic characteristics of the melting 
behavior of metallic microwire and nanowire meshes, it 
would be helpful to find a common parameter which is 
independent of geometrical and physical properties of 
the mesh. In order to deduce such a parameter, let us 
consider a simple model of a wire subjected to a con- 
stant current as shown in Figure 2a. By neglecting the 
difference between and T{y\ t f) for simple approxi- 
mation, the following equation can be easily obtained 
from Equation 4: 



J 8A 



(9) 



where T c is the maximum temperature occurring in the 
center of the wire with x = 112. It indicates that j 2 l 2 (pl\)l 
(T c - T( it j)) is independent of geometrical and physical 
properties of the wire. 

Based on the above consideration, the following dimen- 
sionless parameter Z was proposed as figure of merit of the 
mesh: 



(10) 



which indicates the current-carrying ability of the mesh. 
The variation of calculated Z during the melting process 
is shown in Figure 5, which was developed from the nu- 
merical results in Figure 3. Note that the maximum value 
of Z (i.e., Z c ) corresponding to the maximum value of I m 
(i.e., I mC ) characterizes the current-carrying capacity of the 
mesh, at which the mesh equipped with current-controlled 
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Figure 5 Variation of Z with n b in the melting process of 
both meshes. 



current source will melt until open. Reasonably consistent 
behavior of Z in both meshes can be clearly observed in 
Figure 5. The minor difference can be attributed to the dif- 
ferent melting pathways (see Figure 4), which can be re- 
moved by employing much smaller AI for the microwire 
mesh with sacrifice of computational cost. 

Generally, for the same material, T m , p, A, and A are 
dependent on wire size, while S is dependent on mesh 
structure. For a given mesh structure with a known S, 
the smaller A results in smaller T m and A but larger p, 
and therefore smaller I m according to Equation 10. This 
point is the same with the above numerical results where 
the I m of the microwire mesh is significantly higher than 
that of the nanowire mesh (see Figure 3a). 
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— »— Ag nanowire mesh (Simulation) 
- O- Ag nanowire mesh (Conversion) 
Al nanowire mesh (Conversion) - 



2 4 6 8 10 12 14 16 18 
Melting voltage of wire mesh V m (V) 

Figure 6 Predicted melting behavior of Ag and Al nanowire 
meshes by conversion. 
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Therefore, it is expected that the obtained melting 
behavior of the microwire mesh can be used to predict 
that of the wire mesh with same structure at the same 
working condition even if made from a different wire 
(i.e., different size, different material) through simple 
conversion with the known Z. Taking the Ag nanowire 
mesh as an example, the conversion process is summa- 
rized here. First, the melting current I m for the nano- 
wire mesh can be calculated from Equation 10 with the 
known Z. Second, the variation of the R m for nanowire 
mesh can be calculated from that for the microwire 
mesh in Figure 3b as 



NW 



Rn 



MW x 



NW x 



|MW? 



(ii) 



because of the same melting process. Note that '|nw 
and 'Imw' indicate the case for the Ag nanowire mesh 
and Ag microwire mesh, respectively. Third, the vari- 
ation of V m for the Ag nanowire mesh can be calculated 
by multiplying the obtained R m and I m from the above 
two steps. The predicted melting behavior of the Ag 
nanowire mesh derived from the above indirect conver- 
sion is shown in Figure 6, which indicates good agreement 
with that obtained from direct numerical simulation, and 
therefore validates the feasibility of the present conversion 
method. Figure 6 also gives the predicted melting behavior 
of the Al nanowire mesh with the same structure through 
indirect conversion. Obviously, the melting behavior of 
the mesh is largely dependent on the physical properties 
of the wire itself. 

It should be noted that the present boundary condi- 
tions and mesh structure are only one example. Cer- 
tainly, boundary conditions and mesh structure will have 
great effect on the melting behavior of the wire mesh as 
well as physical properties of the wire itself. However, 
the consistent feature in the melting behavior among the 
wire meshes with the same structure under the same 
boundary conditions will not change. Therefore, the 
present findings can provide meaningful insight for the 
experimental investigation on the reliability of the metal- 
lic nanowire mesh-based TCE. 

Conclusions 

In this work, the melting behavior of an Ag microwire 
mesh induced by Joule heating was numerically inves- 
tigated and compared with that of the corresponding 
Ag nanowire mesh with the same mesh structure but 
different geometrical and physical properties of the 
wire itself. The repetitive zigzag pattern in the rela- 
tionship of melting current and melting voltage dur- 
ing the melting process in the Ag microwire mesh 
was found to be similar with that of the Ag nanowire 
mesh. A dimensionless parameter Z was proposed as 



figure of merit to characterize the current-carrying 
ability of the mesh. The consistent behavior of figure of 
merit in both meshes indicates that the known Z and 
the melting behavior of the Ag microwire mesh can be 
used to predict the melting behavior of the nanowire 
mesh even with different materials (e.g., Ag nanowire 
mesh, Al nanowire mesh), which is hindered by the cost 
of sample preparation and the difficult control of ultra- 
low current stressing in experiments. The present find- 
ings indicate great insight for reliability analysis on the 
metallic nanowire mesh-based TCE, which will be 
beneficial to improve the performance of the corre- 
sponding optoelectronic devices. 
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